Transition from 3D to ID in Bose Gases at Zero Temperature 
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We investigate the effects of dimensional reduction in Bose gases induced by a strong harmonic 
confinement in the transverse cyhndric radial direction. By using a generalized Lieb-Liniger theory, 
based on a variational treatment of the transverse width of the Bose gas, we analyze the transition 
from a 3D Bose-Einstein condensate to the ID Tonks- Girardeau gas. The sound velocity and the 
frequency of the lowest compressional mode give a clear signature of the regime involved. We study 
also the case of negative scattering length deriving the phase diagram of the Bose gas (uniform, 
single soliton, multi soliton and collapsed) in toroidal confinement. 
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I. INTRODUCTION 



Recently the ID regime has been experimentally 
achieved for dilute Bose gases in optical traps at ultra-low 
temperatures by imposing a strong cylindric radial con- 
finement, which inhibits transverse excitations [1]. While 
a 3D dilute Bose gas is well described by the Gross- 
Pitaevskii (GP) theory [2], in the ID regime quantum 
fluctuations become important and an appropriate theo- 
retical description requires the many-body approach in- 
troduced by Lieb and Liniger many years ago [3] and re- 
cently extended by Lieb, Seiringer and Yngvason [4]. In 
particular, at sufRcicntly low densities the Bose gas enters 
in the so-called Tonks-Girardeau (TG) regime, where the 
system behaves as a gas of spinless Fermions (impenetra- 
ble Bosons) [5]. This TG regime has not been achieved 
yet, but is one of the main focus areas of experimental 
investigation. 

In this article we develop a generalized Lieb-Liniger 
(GLL) theory to investigate the crossover from a 3D 
Bose-Einstein condensate (BEG) to a ID TG gas. In 
GLL theory the A'^-body wave function is factorized 
in a cylindric-longitudinal axial wave function and a 
cylindric-radial transverse wave function as recently done 
by Das, Girardeau and Wright (DGW) [6]. While DGW 
have numerically calculated the transverse wave function, 
in our GLL approach the transverse wave function has a 
Gaussian shape with its width as variational parameter. 
In the case of an axially uniform Bose gas the GLL theory 
gives a ground-state energy very close to the one calcu- 
lated in [6]. We use the GLL theory, which is compu- 
tationally much less demanding than DGW method, to 
study dynamical properties of nonuniform gases extend- 
ing our previous analysis of the 3D— >1D transition with 
BEGs [7] by taking into account ID quantum correlations 
and also the effect of an attractive interaction. 



II. BOSE GAS UNDER TRANSVERSE 
CONFINEMENT 

The Hamiltonian of a gas of N interacting identical 
Bose atoms confined in the transverse cylindric radial 
direction with a harmonic potential of frequency ll!±_ is 
given by 
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where V{Yi,Yj) is the inter-atomic interaction. In the 
Hamiltonian we have used scaled units: energies are in 
units of the energy fiuj^ of the transverse confinement 
and lengths in units of the characteristic harmonic length 
a±_ = ^h/{mLLi±). We set V{ri,rj) = gao S^^\ri -rj), 
where g^u — 47ras/a_L with the s-wave scattering 
length of the inter-atomic potential. This Fermi pseudo- 
potential gives the correct dilute gas limit and a well- 
posed variational problem by choosing a smooth trial 
wave function [6]. Note that a confinement-induced res- 
onance (GIR) has been predicted at as/a± ~ 1 [8], i.e. 
when the absolute value \Eb\ of the highest bound-state 
energy Eb — —h^/ima^) of a realistic inter-atomic po- 
tential approaches the confining transverse energy hujj_ = 
h'^/{ma\) [9]. Therefore our analysis will be limited to 
the range |as|/aj^ ^ 10~^, a regime where CIR effects 
are negligible [10]. 

The determination of the A^-body wave function 
'(/'(ri, Tat) that minimizes the energy of the 

system is a very difficult task. Nevertheless, due to the 
symmetry of the problem, a variational trial wave func- 
tion can be written in the form 
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V'(ri, ...,rAr) = /(zi,...,2Ar) Jl 



exp(-^j 



(2) 



where cr is a variational parameter that gives the effective 
transverse length of the Bose gas. For repulsive inter- 
atomic interaction one has a > 1 and only in the case of 
very strong transverse confinement o" = 1 (ID system). 



1 



By rising Eq. (2) and integrating over x and 
y, the total energy per unit of length reads £ = 
£z + where the longitudinal axial energy is £z = 
(/I Ef=i 9iD{cr) Ef<,=i - z^m/L , with 

gioicr) = gsD / C^TTa"^) the effective one-dimensional inter- 
atomic strength, and the transverse radial energy is 
£±_ ~ + cr^)/0, with p — N/L is the homogeneous 

linear density of the Bose gas. 

Due to the Lieb-Liniger theorem [3], for Og > the 
axial energy per unit length £z can be written in terms 
of the axial density p only and the total energy £ can 
thus be rewritten as 



(3) 



where 7 = 2as/a± is the inter-atomic strength and e{x) is 
the Lieb-Liniger function, which is defined as the solution 
of a Predholm equation and it is such that e{x) = x — 
4/(37r)x3/2 for a; < 1 and e{x) = {ir'^ /3){x/{x + 2)f for 
x» 1 [3,11]. 
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Fig. 1. Axial energy per particle Ez/N versus axial den- 
sity p for a Bosc gas under transverse harmonic confinement 
of frequency uj±_ with 7 = 2as/ax = 0.1. Comparison be- 
tween Lieb-Liniger (LL), generalized Lieb-Liniger (GLL) and 
Das-Girardeau- Wright (DGW) theories. The points of DGW 
theory are taken from Fig. 1 of paper [6]. Energy in units 
'hu)^_, length in units ax, and axial density in units 1/ax. 

The energy £ depends on a variational parameter: the 
transverse width a. The minimization of the energy £ 
with respect to a gives the equation 



C7^ = l+7pe'(-^). 



pa^ 



(4) 



Prom this nontrivial equation, for a fixed value of 7, one 
numerically finds the transverse width a of the Bosc gas 
as a function of the longitudinal density p. Analytical 
results can be obtained in limiting cases. In particular, 
under the condition 7 < 1, one gets a = yj^fp for p^ 1/^ 
(3D regime) and cr = 1 for p <C I/7 (ID regime). Note 
that the TG regime is reached at low densities: p 4C 7. 



In Figure 1 we plot the axial energy per particle E^/N 
as a function of the axial density p obtained by using 
different theories: the present generalized Lieb-Liniger 
(GLL), Lieb-Liniger (LL) [3], and Das-Girardeau- Wright 
(DGW) [6]. The results clearly show that GLL approach 
gives practically the same results of the the DGW the- 
ory, that is based on a fully numerical treatment of the 
transverse wave function [6]. Moreover, Figure 1 shows 
that the LL theory, obtained from Eq. (4) simply setting 
cr = 1, is not reliable at high axial densities because the 
Bose gas becomes three-dimensional. 
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Fig. 2. First sound axial velocity ci versus axial density 
p for a Bose gas under transverse harmonic confinement of 
frequency a»x with 7 = 2os/ox = 0.1. Length in units ox, 
time in units wj^ and axial density in units l/ox. 

Having derived the function cr{p), one can numerically 
determine, for a fixed 7, the chemical potential M = f| 
as a function of the density p. It is given by 



2P — ) + Z2 • 



pa- 



a- 



(5) 



Again, analytical results are obtained in limiting cases. 
In particular, under the condition 7 < 1, one finds p = 
\/2(7p)^/^ for p 3> 1/7 (3D regime), p = 7/9-1-1 for 
7 <C (0 < 1/7 (ID regime), and p = Zir'^p^/Q + 1 for 
p <^"1 (TG regime). 

The chemical potential p can be used to derive the ax- 
ial sound velocity ci of the Bose gas, that is obtained 

with the formula ci = \Jp^ [13]- In Figure 2 we plot 
the sound velocity c\ as a function of the density p for 
different values of the inter-atomic strength 7. The re- 
sults show that our generalized Lieb-Liniger (GLL) the- 
ory works remarkably well. It gives the 3D BEG result 
of GP theory at high densities and the TG result of LL 
theory at low densities. 
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III. INCLUDING A LONGITUDINAL 
EXTERNAL FIELD 

Now wc consider the effect of a longitudinal confine- 
ment due to an external potential U{z). Assuming that 
the density p{z, t) varies sufficiently slowly, i.e. the local 
density approximation (LDA) , at each point z the gas can 
be considered in local equilibrium, and the local chemical 
potential n[p{z, t)] is provided by Eq. (5) supplemented 
by Eq. (4). Under LDA the dynamics can be described 
by means of a complex classical field ^{z, t) that satisfies 
the time-dependent GLL nonlinear Schrodinger equation 

where p{z,t) = This equation is a Euler- 

Lagrange equation obtained by imposing the least-action 
principle to the following GLL action functional 

(7) 

where £ is given by Eq. (3), and minimizing with respect 
to The other Euler-Lagrangc equation is precisely 
Eq. (4), that can be obtained from A by minimizing 
with respect to a. It is important to observe that at 
high densities (p ^ I/7) A gives the action functional 
of the non-polynomial Schrodinger equation (NPSE) we 
have derived to describe BECs under transverse confine- 
ment [7], while at low densities (p -C I/7) and in the 
static case A reduces exactly to the energy functional 
recently introduced by Lieb, Seiringer and Yngvason [4]. 
Note that Eq. (6) is equivalent to the two hydrodynamics 
equations of a viscousless fluid in ID with density field 
p{z,t) and velocity field v{z,t) = (ggg for 

instance [12]). 

Neglecting the kinetic gradient term |^jt (Thomas- 
Fermi approximation) and imposing the equilibrium con- 
dition, we obtain the equation 

p[po{z)] =iiT- U{z) , (8) 

where pq{z) is the equilibrium density and pT is the 
chemical potential of the axially inhomogeneous system. 
Inverting Eq. (8) one obtains po{z) = p~^[pT — U{z)]. 
The value of px is found by imposing the normalization 
condition N = J pQ{z)dz. In the case of axial harmonic 
confinement with scaled potential U{z) = \}?z'^, where 
A = lOz/lo±^ is the anisotropy of the full harmonic trap 
with Uz axial harmonic frequency, the normalization con- 
dition can be rewritten in the form 

where Q = zX/^/2pT is the rescaled axial coordinate. 
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Fig. 3. Scaled axial density profile po(C)/Po(0) versus 
scaled axial coordinate (" = z\/^2pT for a Bose gas of A'^ 
atoms in a harmonic confinement with 7 = 2as/a± = 0.1. 
A = u>z/ijJ± is the trap anisotropy and pr is the chemical 
potential. Units as in Fig. 1. 

In Figure 3 we plot the axial density poiO- oi'" 
der to compare different cases, we rescale the density 
with respect to the central one po{0)- We have verified 

that po{C) follows with good accuracy the law po{C) = 
po{0){l — C^)") where the exponent a ranges between 
a = 2 (TG regime) to a = 1/2 (3D BEG regime) by 
increasing A'^A. 

The collective oscillations arc determined by writing 
the density field in the form p{z, t) — po{z) + e"^^^ pi{z), 
with the function pi{z) obeying a linearized equation 
which follows from hydrodynamic equations [12,13]. 
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Fig. 4. Frequency CI of the lowest compressional mode of a 

Bose gas of N atoms in a harmonic confinement. A = lOzI^i_ 
is the trap anisotropy and 7 = 2as/a± is the inter-atomic 
strength. Units as in Fig. 2. 

The problem has simple analytic solutions if the equi- 
librium density po{z) has the form po{z) = po{0){l — 
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z"^ /Z"^)"-. In this case the lowest axial compressional 
mode is given by 



For a fixed value of 7 we easily determine a as a func- 
tion of iVA using Eq. (4), Eq. (5), Eq. (8) and Eq. (9). 
In such a way we obtain the lowest axial compressional 
mode il as a function of NX, as shown in Figure 4. As 
expected, while for a large 7 there arc two plateaus cor- 
responding to ri^/w^ = 5/2 at high densities (3D GP) 
and n^/wf = 4 at low densities (TG), for A < 1 an ad- 
ditional plateau appears with O^/w^ = 3 at intermediate 
densities (ID GP). 

The results shown in Figure 4 improve those obtained 
by Menotti and Stringari [12] using a sum-rule approach 
in both LL theory (ID regime) and GP theory (3D BEG 
regime). In Ref. [12] they found that the curves of f^^/w^ 
obtained with the two theories do not match at interme- 
diate densities if 7 ~ 0.1. In fact, LL theory predicts the 
achievement of the ID mean-field regime when GP theory 
already exhibits significant 3D effects. As clearly shown 
in Figure 4, GLL theory overcomes such difficulties. 

Here we have shown that the time-dependent LDA re- 
produces known results in the appropriate limits and 
therefore can be safely used to calculate collective os- 
cillations (one-phonon excitations) around static config- 
urations of the Bosc gas. Note that, in general, time- 
dependent LDA can give incorrect results for more com- 
plex dynamical phenomena, like multi-phonon excita- 
tions [14] or the interference of two Bosonic clouds [15]. 

IV. NEGATIVE SCATTERING LENGTH 

Our variational approach can be applied also to the 
case of negative 3D scattering length (a^ < 0) by using 
Eq. (6) and Eq. (7) and simply setting e{x) = x with 
a; < in Eq. (3), Eq. (4) and Eq. (5). This corresponds 
to the Hartree approximation HjLi ^(-^j) of the the ex- 
act N-body wave function /(zi, zjv) [7]. In the pure 
ID regime (cr = 1), Calogero and Degasperis [16] have 
shown that with negative scattering length the exact N- 
body wave function is well approximated by the Hartree 
ansatz. In this way Eq. (6) reduces to the NPSE, Eq. 
(4) gives the equation = y/1 — \'y\p for the transverse 
width [17]. 

In a toroidal configuration, i.e. cylindric confinement 
with periodic boundary conditions in the axial direction 
and U{z) = 0, the Bose gas with negative scattering 
length is axially uniform below a critical number Nc of 
Bosons and becomes localized (bright soliton) only above 
this critical number. Garr, Glark and Reinhardt [18] have 
shown that in the pure ID regime {a = 1) this critical 
number is given by Nc = tt^/ (|7|L), where L is the length 
of the torus. We easily extend that prediction to a 3D 
Bose gas by calculating with NPSE the wave-length Ac of 



the zero-energy Bogoliubov excitation above the uniform 
state [19]: in this way, setting Ac — L, we find the cqua- 
tion Ar,|7|L/7r2 = (1 - iV,|7|/i)-V2/(i _ 3A^c|7l/(4i))- 
Gontrary to the ID theory, our GLL theory predicts the 
existence of the single bright soliton in toroidal confine- 
ment only below a critical strength A''|7|. 
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Fig. 5. Phase diagram of a Bosc gas with negative scat- 
tering length in toroidal confinement. N is the number of 
Bosons, 7 = 2as/a± is the inter-atomic strength and L is the 
length of the torus. Units as in Fig. 1. 

We have numerically verified that such a critical 
strength is well represented by A''|7| = 4/3. Above this 
line only multi-soliton configurations are possible, which 
collapse at A^|7| = L/2 as in the case of the NPSE with 
box boundary conditions (see [17]). In Figure 5 we show 
the phase diagram of the toroidal Bose gas with negative 
scattering length in the plane A'^|7| vs L: uniform, local- 
ized with a single bright soliton, localized with a multi 
soliton train and collapsed. The properties of these multi- 
soliton configurations will be discussed elsewhere. 

V. CONCLUSIONS 

We have introduced a very effective generalized Lieb- 
Liniger theory, based on a variational treatment of the 
transverse width, that points out the crucial interplay 
between quantum correlations and dimensionality in con- 
fined Bose systems. The first sound velocity, the density 
profile and the collective oscillations give a clear signature 
of the regime involved. Our results exactly reproduces 
the ones obtained in specific regimes with the hydrody- 
namic method or the sum-rule method, but, they also 
give the crossover from one regime to another. Finally, 
we have applied our approach to a Bose gas with negative 
scattering length deriving the complex phase diagram 
(that includes uniform, single-soliton and multi-soliton 
configurations) of an attractive Bose gas in toroidal con- 
finement. 
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